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ABSTRACT 

A major question in the study of black hole binaries involves our understanding of the accretion geometry 
when the sources are in the "hard" state. In this state, the X-ray energy spectrum is dominated by a hard power- 
law component and radio observations indicate the presence of a steady and powerful "compact" jet. Although 
the common hard state picture is that the accretion disk is truncated, perhaps at hundreds of gravitational radii 
(/?g) from the black hole, recent results for the recurrent transient GX 339-4 by Miller and co-workers show 
evidence for optically thick material very close to the black hole's innermost stable circular orbit. That work 
focused on an observation of GX 339^ at a luminosity of about 5% of the Eddington limit (Z-Edd) and used 
parameters from a relativistic reflection model and the presence of a soft, thermal component as diagnostics. In 
this work, we use similar diagnostics, but extend the study to lower luminosities (2.3% and 0.8% Lndd) using 
Swift and RXTE observations of GX 339^. We detect a thermal component with an inner disk temperature 
of ^0.2 keV at 2.3% Lndd- At 0.8% Lndd, the spectrum is consistent with the presence of such a component, 
but the component is not required with high confidence. At both luminosities, we detect broad features due 
to iron Ka that are likely related to reflection of hard X-rays off the optically thick material. If these features 
are broadened by relativistic effects, they indicate that optically thick material resides within 10 Rg down to 
0.8% Lfidd, and the measurements are consistent with the inner radius of the disk remaining at ^4 Rg down to 
this level. However, we also discuss an alternative model for the broadening, and we note that the evolution 
of the thermal component is not entirely consistent with the constant inner radius interpretation. Finally, we 
discuss the results in terms of recent theoretical work by Liu and co-workers on the possibility that material 
may condense out of an Advection-Dominated Accretion Flow to maintain an inner optically thick disk. 

Subject headings: accretion, accretion disks — black hole physics — stars; individual (GX 339-4) — X-rays: 
stars — X-rays: general 



1. INTRODUCTION 

Over the past several years, progress has been made in con- 
straining the geometry of black hole accretion disks. Most 
binaries with stellar mass black holes are X-ray transients, 
and X-ray observations made when the sources are bright 
(Lx ^ 10^^"^'' ergs s"') have uncovered iron Ka emission lines 
with broad and redshifted profiles thought to be produced via 
fluoresce nce when hard X-rays refle ct off optically thick disk 
material ([Reynolds & NowaHl2003l) . Fitting the energy spec- 
tra with relativistic reflection models yield inner disk radii 
close to the innermost stable circular orbit (ISCO) of the black 
hole ( Miller et al. 2002). Another constraint comes from 
the high frequency (100-500 Hz) quasi-periodic oscillations 

i QPOs) that are prese nt for a number of black hole systems 
Remillard et al.ll2002h . Although the origin of these QPOs is 
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Still unclear, if they are related to Keplerian orbital frequen- 
cies f or matter in an optically thick (e.g., Sha kura & SunvaevI 
Il973h accretion disk, they correspond to time scales very 
close to the ISCO, implying that the disk must extend near 
the black hole. Furthermore, these constraints on the inner 
radius of the accretion disk (Rm) occurred at times when the 
X-ray energy spectrum in cluded a strong therrnal co mponent 
that is consistent with the IShakura & SunvaevI d 1973b model 
with an inner disk temperature of ~ 1 ke V and an inner radius 
close to or at the ISCO. 

As black hole outbursts evolve, they enter different spectral 
states, and one definition for the various sta tes was recently 
described in lMcCUntock & Remilla rd (2006). While the iron 
lines, QPOs, and thermal components discussed above have 
provided a determination of Ri„ in the "thermal-dominant" 
and "steep power-law" (SPL) spectral states, in this work, 
we focus on the hard state where the accretion disk geom- 
etry is still unclear. The hard state is most often observed 
at the beginning and end of the outburst and is characterized 
by a hard energy spectrum and a high level of X-ray timing 
noise. At the ends of outbursts, the hard state is most often 
seen when sources reach Eddington-scal ed luminosities be- 
low L/LEdd = 0.01-0.04 (Maccaron ell2003l) . but the hard state 
also can occur at higher luminosities, especially when out- 
bursts commence. Radio observations have made it clear that 
a steady outflow in the form of a "co mpact" jet is character- 
istic of the hard state (lFendeiir2001h . and, at least in some 
systems, the jet power is inferred to be more than 20% of the 
X-ray luminosity (Fender et al. 2001). Thus, one reason that it 
is important to understand the disk geometry and other physi- 
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cal processes that occur in the hard state is that this is the one 
state that is linked to the production of a steady and powerful 
jet. 

The common picture for the hard state accretion 
disk has been that its inner edge recedes (i.e., Rin 
increases), leaving a hot flow, such as an Advection- 
Dominated Accretion Flow (ADAF) or a spherical 
corona, where most of the X-ray emissi on is p r oduced 
via Comptonization i Naravan. McClinto ck & Yil 119961 : 
lEsin. McCUntock & N aravan J997l) . Indeed, it has 
been shown that evaporation of the accretion disk can 
lead to a truncated disk at low m ass accretion rates 
(iMeyer. Liu & Meyer-Hofmeisted l2000l) . From an obser- 
vational standpoint, several measurements are suggestive 
of a truncated disk, but no one method provides a direct 
measurement of Ri„. Examples of measurements that support 
an increase in include a rapid drop in the temperature 
and flux of the thermal component (often to non-detection in 
the >3 keV band) as the source enters the hard state as well 
as a drop in the characteristic frequencies seen in the power 
spectra ( Poutanen, Krolik & Rvde 1997; Doveetal. 1997; 
Tomsick & Kaaret 2000; Revnivtsev, Gilfanov & Churazqy 
200 k .Kalemci et aL ,2004; Tomsick. Kalemci & Kaarea 
20041) . In addition, a gradual drop in the strength of the 
reflection component can be explained if the disk is trun- 
cated and there is over lap betwe en the disk and cor ona 
jZdziarski. Lubinski & Sm ith 1999; IZdziarski et al.ll2003l) . 

While these measurements can be explained within the 
truncated disk picture, other explanations have also been sug- 
gested. Much of the evidence for the disappearance of the 
thermal component come from RXTE spectra that lack soft 
X-ray coverage, allowing for the possibility that the thermal 
component simply shifts out of the instrumental bandpass. 
Also, an increase in the ionization state of the inner disk 
could cause a drop in the strength of the reflection compo- 
nent. In addition, some recent X-ray observations of black 
holes in their hard state cast doubt on whether the truncated 
disk picture is correct. The black hole system GX 339^ 
was observed in the hard state with XMM-Newton and the 
Rossi X-ray Timing Explorer (RXTE), and a fit with reflec- 
tion and iron line emission components accounting for rela- 
tivistic smearing effects are consist ent wit h an inner disk ra- 
dius very close to the ISCO (Miller et alJl2006) . The ISCO 
is at 6^g (where Rg = GMjc^, G and c are constants, and M 
is the black hole mass) for a non-rotating black hole and at 
1.23/?g for a black hole ro tating at the limiting spin rate of 
fl* = 0.9982 (T horneiri974l) . w here a^, is the mass-normalized 
angular momentum parameter. iMiller et"an (12006 ) obtained a 
value of /?i„/«g = 4.0 ± 0.5 from the GX 339-4 hard state X- 
ray spectra. As this value is only s lightly larger than the value 
of 2-3 /?g that MillefetaD dlOOl obtained for GX 339-4 in 
the SPL state, it suggests that there is very little change in 
7?in between the two states. However, we must keep in mind 
that the interpretation of the relativistic smearing parameters 
depend on the model being physically correct, and other ex- 
planations have been adv anced for iron line broadening (e.g., 
iLaurent & T itarchuk 2007 ) . 

For the thermal component from an optically thick disk, one 
expects the emission to fall outside the X-ray regime if Ri^ in- 
creases dramatically (e.g., to hundr eds of Rc , or more) as pre- 
dicted by ADAF models (McCHnto ck et al.l l2b03). However, 
there are numerous examples of black hole systems, includ- 
ing well-known systems such as Cygnus X-1 and GX 339^, 
that exhibit thermal X-ray emission with temperatures 0.1- 



0.4 keV when they ar e in the hard state (lEbisawa et al.|[T996t 
IZdziarski et al.l Il998l) . While this implies that very large 
values of Ri„ are not required for black holes to enter the 
hard state, the thermal components have been found mostly 
only when the systems are in the brightest phases of their 
hard states. For example, the XMM-NewtonlRXTE obser- 
vation of GX 339-4 in the hard state described above oc- 
curred at the start of an outburst at L/L-em = 0.05 and spec- 
tral fits yielded a thermal component with a temperature of 
^^0.4 keV. X-ray observations of other black hole systems at 
lower X-ray luminosities have provided relatively high qual- 
ity X-ray spectra, a nd, in some cases, soft compone nts have 
been detected (e.g., iMiller. Homan & Miniutti"2006'). while, 
in other cases, they have not (e.^ ., McClintock et al. 20011 
iTomsick. Kalemci & Kaar et 2004). 

For this paper, we have obtained broadband Swift 
(Gehrels et al. 2004) and RXTE (iBradt. Rothschild & SwaiS 
1 19931) observations of the recurrent transient GX 339^ in 
the hard state to constrain the accretion disk geometry in this 
state. The observations were made at the end of th e most re- 
cent outburst that began in 2007 January (iMiUer et al. 2007|) 
after the transition to the hard state that occurred in 2007 May 
(iKalemci et al.ll2007h . The luminosities of the observations 
we use in this work are well below L/LEdd = 0.05, and our goal 
is to use the spectra along with reflection models that account 
for relativistic effects to constrain the accretion geometry at 
these low luminosities. 

2. OBSERVATIONS 

2.1. RXTE Monitoring Observations 

We obtained daily pointed RXTE monitoring observations 
of GX 339-4 starting on 2007 April 20 (MJD 54,210) when 
the flux dropped below 4 x 10"^ ergs cm"- s"' as measured 
in the 1.5-12 keV band by the RXTE All-Sky Monitor. Ob- 
servations were made under our program (proposal #92704) 
until 2007 July 18 (MJD 5 4,299), and the RX TE Propor- 
tional Counter AiTay (PCA iJahoda et all l2006h 3-25 keV 
light curve during this 89 day period is shown in Figure [T] 
The exposure times were typically 1-3 ks for these moni- 
toring observations although some longer RXTE observations 
were also obtained as described below. We reported an in- 
crease in timing noise and a hardening of the GX 339-4 en- 
ergy spec trum that indicates a chan ge to the hard-intermediate 
state (seelHonian & Bellonill2005h occurred on 2007 May 12 
(IKalemci etal.ll2007h . and this is marked in Figure [T] with a 
vertical dotted line at MJD 54,232. By 2007 May 22 (marked 
in Figure [T] by the vertical dashed line at MJD 54,242), the 
source had reached the hard state with an energy spectrum 
dominated by a power-law with a photon index of F = 1.6. 
Further evidence that the source reached the hard state by 
this time inc ludes an increase in the optical and infrared flux 
(iBuxton & B ailyn 2007) as well as core radio flux (S. Corbel, 
private communication), which is an indication for the pres- 
ence of a compact jet. 

2.2. Swift one/ RXTE Observations 

The main focus of this work is the study of the broadband 
~ 1-200 ke V energy spectra from GX 339^ at two times after 
the source made a transition to the hard state. As indicated in 
Figure[T] the first Swift observation occurred on 2007 May 25 
(MJD 54,245), a few days after the transition to the hard state. 
The average Swift X-ray Telescope (XRT Burrows et al...2005i) 
count rate during the 6150 s observation was 10.41 ±0.04 
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Fig. 1. — The 3-25 keV light curve for GX 339^ as measured by the 
RXTE/PCA. The rates shown are for a single Proportional Counter Unit 
(PCU2). The vertical dotted line shows the start of the transition to the hard 
state, and the vertical dashed line shows when the source reached the hard 
state. The times of the Swift observations at MJD 54,245 and MJD 54,261- 
54,265 are marked. 

c s"' (0.8-8 keV). As indicated in Table[T] we obtained RXTE 
observations that were simultaneous with much of the Swift 
observation. The average PCA count rate was 44. 1 ± 0. 1 c s"' 
(3.6-25 keV) for Proportional Counter Unit (PCU) #2, and 
the average High Ener gy X-ray Timing Experiment (HEXTE 
iRothschild e"tal]ll998h count rate was 10.4 ±0.2 c s"' (17- 
240 keV) for HEXTE cluster B. As described below, Spec- 
trum #1 consists of the spectra from XRT, PCA, and HEXTE 
from these observations. 

We used the same 3 instruments for Spectrum #2, but the 
observations were made 2-3 weeks later during the time pe- 
riod 2007 June 10-14. Table[T]shows that there were 3 Swift 
and 5 RXTE observations made during this time frame. Al- 
though the Swift and RXTE observations were mostly not 
strictly simultaneous, the source did not show large changes in 
count rate or spectral hardness over the 5 day period. The 3.6- 
25 keV PCA count rates for each of the 5 observations listed 
in Table [U were 16.4±0.1 c s"', 16.2±0.1 c s"', 15.7 ±0.2 
c s-i, 15.1 ±0.1 c s"', and 14.2±0.1 c s"', showing a grad- 
ual but not dramatic decline. As shown in Figure [T] this pe- 
riod is close to the minimum flux that GX 339^ obtained be- 
fore its flux began increasing again. The average Swift XRT 
count rate was 2.86 ±0.02 c s"' (0.8-8 keV), while the av- 
erage PCU2 and HEXTE-B count rates were 15.6 ±0.1 c s"' 
(3.6-25 keV) and 3.8 ±0.2 c s"' (17-240 keV), respectively. 
Thus, when compared to Spectrum #1, the XRT count rate is 
lower by a factor of 3.6 while the count rates for the RXTE 
instruments are lower by a factor of ^2.8. 

3. SPECTRAL ANALYSIS 

We performed spectral analysis using XRT, PCA, and 
HEXTE data from the observations described above. To ex- 
tract the spectra, we used the Swift and RXTE tools provided 
in the HEASOFT v6.3.1 software package. The XRT instru- 
ment consists of a CCD imager at the focus of a grazing in- 
cidence X-ray telescope. Although the instrument is capable 
of two-dimensional imaging for faint sources, our target was 
bright enough to require the use of Windowed Timing mode, 
which provides one-dimensional imaging. For the GX 339- 
4 energy spectra, we determined the one-dimensional cen- 
troid of the target and extracted the photons within 47" of 
the centroid. We also extracted background spectra that in- 



clude photons from regions between 77" and 171" from the 
GX 339^ centroid. We used version 9 of the XRT re- 
sponse matrix and kept events with grades in the range 0- 
2 (swxwt0_20010101v00 9 . rmf). In addition, we pro- 
duced an ancillary response file (ARE) with the HEASOFT 
tool xrtmkarf vO.5.3. In making the ARE file, we 
used an exposure map produced with the tool xrtexpomap 
vO.2.2. 

For the PCA spectra, we used the SkyVLE background 
model that was most recently updated 2005 November 28^, 
and we used the HEASOFT tool pcarsp vlO . 1 to pro- 
duce the response matrix. Although the PCA has 5 PCUs, it 
is typical that only 2 or 3 PCUs will be turned on for a given 
observation. In our case, only PCUs and 2 were turned on 
for all of our RXTE observations. Due to the loss of the PCU 
propane layer, this PCU's response is not as well known as 
the other PCUs (Jahoda et al. 2006). Also, for all the PCUs, 
the response for the top anode layer (where most of the counts 
are detected) is better modeled than the bottom layers. Thus, 
for this study, we only used the top anode layer of PCU 2 
when extracting spectra. To check on the level of systematic 
error, we extracted the PCU 2 spectrum for the Crab nebula 
using Observation ID 92802-01-22-00, which is an observa- 
tion from 2007 May 14. Fitting the spectrum with an absorbed 
power-law, we find that 1 % systematic errors are required to 
reach a reduced-x^ near 1 .0, and we use 1 % systematics for 
the GX 339^ spectra as well. For the HEXTE spectra, we 
used only HEXTE-B because HEXTE- A no longer obtains 
background measurements. 

We performed preliminary spectral fits to Spectrum #1 to 
check on the agreement between the calibrations of the dif- 
ferent instruments. For these preliminary fits, we used XRT, 
PCA, and HEXTE data in the 0.3-10 keV, 2.75-25 keV, and 
1 7-240 keV energy bands, respectively. We used XSPEC vl2 
and fitted the spectra with an absorbed power-law model. We 
also included a multiplicative constant in the model to allow 
for differences in the overall normalizations between instru- 
ments. The fact that we obtain a poor fit (x^ /v = 1282/269) 
is partially due to the fact that the power-law model is too 
simple as well as being due to instrument calibrations. For 
example, it is known that the XRT calibration is complicated 
by featur es due to the Si02 lay er in the CCD detectors at low 
energies (lOsborne et alJl2005l) . In fact, we see a sharp dip in 
the XRT residuals near 0.5-0.6 keV that is likely instrumen- 
tal. Furthermore, we see that the XRT calibration does not 
agree with the PCA calibration above 8 keV, and these issues 
lead us to use the XRT in the 0.8-8 keV energy range for the 
fits described below. In addition, the PCA calibration does 
not agree with XRT at the very bottom of the PCA range as 
there are strong positive residuals in the PCA spectrum be- 
low 3.6 keV. Thus, in the following, we use the 3.6-25 keV 
PCA spectrum. PCA and HEXTE match well in the region 
where they overlap, and we use the full 17-240 keV HEXTE 
bandpass. In addition, we note that there is good agreement 
(within 5%) between the overall normalizations of the three 
instruments. 

For PCA, Spectra #1 and #2 have high statistical qual- 
ity (hundreds or thousands of counts per energy bin) across 
the bandpass, but some rebinning was required for XRT 

' A problem with the software that produces PCA background spectra was 
recently announced (see http://heasarc.gsfc.nasa.gov/docs/xte/xte_lst.html I. 
We re-extracted the spectra with the connected software and found that the 
problem did not impact our observations. 
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and HEXTE. For XRT, we rebinned from 719 channels to 
148 channels, leaving averages of 432 counts/bin and 259 
counts/bin for Spectra #1 and #2, respectively. We rebinned 
the HEXTE spectra from 210 channels to 22 channels. 

4. RESULTS 

4.1. Fits with Basic Models: Is a Thermal Disk Component 

Present? 



We fitted the spectra with 
law model. To account for 
the photoelectric absorption 
iBalucinska-Church & McCam monI 



in absorbed power- 
absorption, we use 
cross sections from 

^ . (119921) and e lemental 

abundances from [Wilms, Allen &^cCray| (|2000|) . which 
correspond to the estimated abundances for the interstellar 
medium. A power-law model does not provide a good fit 
to either spectrum, with equal to 651.8 and 321.6 for 
Spectra #1 and #2, respectively (for 213 degrees of freedom 
in both cases). Figures |2] and |3] show the spectra and the 
residuals in the form of a data-to-model ratio. In both cases, 
the strongest feature in the residuals occurs in the region near 
the iron Ka region. Positive residuals are present close to the 
^6-7 keV range where an emission line might be expected 
to be present, and negative residuals are seen from 7 keV to 
beyond 10 keV. The positive residuals that are present in the 
20-40 keV range, especially in Spectrum #1, may be due to a 
reflection component. Finally, the curvature in the residuals 
at the lower end of Spectrum #1 may indicate the presence of 
a thermal component from an optically thick accretion disk. 
The fit parameters are given in Table |2] and it is notable that 
column density of A^h = (3.1 ±0.1) x 10'' cm"^ obtained for 
the power-law fit to Spectrum #1 is somewhat lower than the 
value inferred from t he work of iHynes et al.l (l2004l) . From 
optical observations, iHvnes et"aL (|2004|) prefer a value of 
E(B-V) ^ 0.85, which corresponds to A^h Z 4.7 x 10^' 
cm~^ using conversions given in iPredehl & Sciimiti (fT995h. 
Thus, the lower column density could also be an indication 
that a thermal component is present. 

To test whether a thermal component is present in the spec- 
tra, we performed the fits detailed in Table |2] For Spectrum 
#1, when we add a disk-blackbody (Mitsuda et al. 1984) com- 
ponent to the power-law, the quality of the fit shows a large 
improvement to x^/i^ = 532.0/211. This, in addition to the 
fact that the A^h increases to a level which is consistent with 
the lower limit from optical extinction measurements, are in- 
dications for the presence of a thermal component. However, 
even with the disk-blackbody component, the fit is poor due to 
the iron features in the spectrum, so we re-fitted Spectrum#l 
with a model that takes these features into account. Although 
we use a more physical model for the iron line and reflec- 
tion below, here we add a smeared iron edge (Ebisawa et al. 
IT994i) to the model because it is a simple addition that signifi- 
cantly improves the fit. When the disk-blackbody component 
is added to a model with a power-law and a smeared edge, 
the fit shows a large improvement from x^/ v = 444.7/21 1 to 
338.3/209. 

We carried out the same series of fits for Spectrum #2, and 
while they also provide evidence for the presence of a ther- 
mal component, the evidence is considerably weaker than for 
Spectrum #1. For the models without the smeared edge (see 
Table [2]i, adding the disk-blackbody component gives a rela- 
tively large improvement in the fit from x^/i^ = 321. 6/213 to 
280.1/21 1. However, with the smeared edge, the change from 
X^/i^ = 210.6/21 1 to 200.2/209 is rather small. 
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Fig. 2.— Spectrum #1, including data from Swift/XKT (0.8-8 keV), 
RXTE/PCA (3.6-25 keV), and RXTE/HEXTE (17-240 keV), fitted with an 
absorbed power-law model. The data-to-model ratio is shown in the bottom 
panel. 
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Fig. 3.— Spectrum #2, including data from the XRT, PCA, and HEXTE 
instruments, fitted with an absorbed power-law model. The data-to-model 
ratio is shown in the bottom panel. 

To further test the significance of the thermal component in 
the spectra and to determine if there is any evolution in the 
thermal component between Spectra #1 and #2, we produced 
error contours for the disk-blackbody temperature (kTin) and 
normalization (A^dbb)- In Figure |4] the outermost contour 
for each spectrum corresponds to 90% confidence for two- 
parameters of interest (Ax^ = 4.61). The error region for 
Spectrum #1 shows that A^dbb is significantly different from 
zero, consistent with the presence of a thermal component 
in Spectrum #1. The error region for Spectrum #2 is well- 
separated from that of Spectrum #1 , showing a clear change in 
the thermal component between the two spectra. The results 
indicate a drop in kTin or A^dbb or both parameters. To estimate 
the significance of the thermal component in Spectrum #2, 
we adjusted Ax^ until the confidence contour reached a A^dbb 
value of zero. The 99% confidence contour (Ax^ = 9.21) does 
not reach zero, but zero is reached with a slightly larger con- 
tour (Ax^ = 10), which is consistent with a ^1% chance that 
this component is spurious. 

4.2. Iron Line and Reflection Modeling 
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Fig. 4. — Confidence contours for the thermal (disk-blackbody) compo- 
nents for Spectra #1 (soHd contours) and #2 (dashed contours). In each case, 
the inner-most contour encircles the 68% confidence (Ax' = 2.30) error re- 
gion for the two parameters and the outer-most contour corresponds to 90% 
confidence (Ax^ = 4.61). 

We re-fitted the spectra in an attempt to improve our mod- 
eling of the iron line and reflection. Initially, we removed the 
smeared edge and fitted the spectra with a model consisting 
of a disk-blackbody, a power-law, and a Gaussian emission 
line. Although the Gaussian greatly improves the fit in both 
cases, the energy is well below the 6.4-7.1 keV iron range 
with values of Eune = keV and 4.0^? j keV for Spectra 
#1 and #2, respectively. As suggested by the residuals in Fig- 
ures |2] and [3] the Gaussians are also very broad with widths 
of a ^ 2.0 keV in both cases. This suggests that we may be 
seeing relativistically smeared iron lines.^ 

The laor emission line model (Laorl l99ll) is appropriate 
for reflection from an accretion disk around a rotating black 
hole, and if the iron line is due to reflection, then one expects 
to see evidence for excess emission from reflection at higher 
energies (^20-40 keV) as well. In fact, the positive residuals 
for Spectrum #1 (and to some extent for Spectrum #2) in this 
energy range suggest that there is indeed a reflection compo- 
nent. Thus, we re-fitted the spectra with a model consisting 
of a disk-blackbody, an iron emission Une, and a p exriv re- 
flection component (Mag dziarz & Zdziarskil[T995h . which in- 
cludes both direct emission from a power-law as well as emis- 
sion reflected off an accretion disk with neutral or partially 
ionized material. While the pexriv model includes the ion- 
ization effects, which are likely important for the hot accre- 
tion disks around X-ray binaries, it does not include relativis- 
tic smearing. Thus, in our model, we convolved the pexriv 
model with the laor model shape using the XSPEC con- 
volution model kdblur. We set up the XSPEC model so 
that kdblur convolves the sum of a narrow iron line and the 
pexriv model. The kdblur parameters include Rin (the in- 
ner radius of the disk), /?out (the outer radius of the disk), / (the 
binary inclination), and q (the power-law index for the radial 
emissivity profile). 

The results of these fits are given in Table |3] and the fit- 
ted spectra are shown in Figures |5] and |6] The power-law and 
disk-blackbody parameters have similar values to those that 
we found with the basic models (Table |2]i, and the quality of 
the fits is better than we obtained with the basic models. For 
both spectra, the presence of the pexriv reflection compo- 



nent is required at high confidence as indicated by the fact that 
the reflection covering fraction (r2/27r) is significantly differ- 
ent from zero. The exact value of ft/lTT depends strongly on 
the binary inclination, which is not known for GX 339^. Pre- 
vious fits to highe r quality X-ra y spectra have given a value of 
i = 20^1^ degrees (iMiller et al.ll2006 ). and we adopt a value of 
20° to facilitate comparisons to the previous work. With this 
inclination, we obtain values of r2/27r = 0.22;|;q [J5 and 0.24;|;q 
for Spectra #1 and #2, respectively. The pexriv ionization 
parameter, ^, is not very well-constrained, but it is signifi- 
cantly greater than zero for both spectra, indicating a disk that 
is at least partially ionized. 

While the reflection component is statistically significant in 
both Spectra #1 and #2, an additional iron line in the 6.4- 
7.1 keV range (see Enn^ in Tabled is required for Spectrum 
#1 but is required at only slightly more than 90% confidence 
for Spectrum #2 as shown by the values of the emission line 
normalization (A^iine) given in Table [3] This is not due to the 
lack of iron features in Spectrum #2 but because, with the 
relativistic broadening, the reflection component contains a 
bump related to the iron absorption edge that can mimic a 
broad iron emission line (see Figures|5]and|6]l. However, even 
though the emission line is not clearly detected in Spectrum 
#2, the parameters that account for the relativistic smearing 
are still well-constrained for both spectra because both the 
line and the reflection (pexriv) components are smeared. 

While there are four relativistic smearing parameters, only 
two of the parameters are left as free parameters in our fits. 
As mentioned above, we fixed the binary inclination to 20°, 
and we fixed the outer disk radius to /?oiit = 400/?g(= GM/c^, 
where G and c are constants and M is the black hole mass). 
One of the free parameters is the inner disk radius, and we 
find that Rin = 3.6t\ f^Rg and 2.9^(^:^7?^ for Specti-a #1 and #2, 
respectively, implying that the reflecting material is very close 
to the black hole. The other free parameter is the power-law 
index for the radial emissivity profile, and we obtain values 
consistent with q = 3 for both spectra. 

5. DISCUSSION 

5.1. Constraints from the Reflection Model 

Our results on the iron Une and reflection component of 
GX 339^ join a relatively small number of observations 
where reflection models that account for the relativistic ef- 
fects near black holes have been fitted to broadband X-ray 
spectra of black holes in the hard state. Our findings are 
most directly comparable to the results of Miller et al. (2006]) 
where both the GX 339^ and Cygnus X-1 showed evidence 
for a broad iron Ka emission line and a smeared reflection 
component while they were in the hard state. Table [3] com- 
pares the parameters for GX 339-4 from the fits to the 2004 
XMM-NewtonlRXTE specti'um reported in Miller et al. (200^ 
to the parameters we obtain by fitting the same model to our 
Swift/RXTE spectra. Independent of assumptions about the 
distance to the source. Spectra #1 and #2 were taken when 
the GX 339^ luminosity (1-100 keV, unabsorbed) was, re- 
spectively, 2.4 and 7.0 times lower than the 2004 spectrum. 
Adopting a source distance of 8 kpc (jHv nes et al. 2004) and 
a black hole mass of 5. 8 Mrr^ (fHvnes et al. .2003) , which are 
the same values used by lMiller et alJ (l2006h . we estimate that 
the Eddington-scaled luminosities during our observations are 
L/LEdd = 0.023 and 0.008. 

The two most significant reflection parameters for answer- 
ing the question of the accretion geometry for GX 339-4 in 
the hard state are the covering fraction ft /In and the inner 
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0.01 1 , , .1 j. ,1 1 . I 

1 10 100 

Energy (keV) 

Fig. 5. — Spectrum#l fitted with the model detailed in Table[3]and plotted 
in flux units. The various model components are shown and include a thermal 
disk-blackbody component, a power-law, a reflection component, and an iron 
emission line (the last two include relativistic effects). 
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0.01 1 . . .i/ ' . , . f , = 

1 10 100 

Energy (keV) 

Fig. 6. — Spectrum #2 fitted with the model detailed in Table[3]and plotted 
in flux units. The various model components are shown and include a thermal 
disk-blackbody component, a power-law, a reflection component, and an iron 
emission line (the last two include relativistic effects). 

radius of the optically thick disk (ri„ = Rin/Rg). The covering 
fraction of ~0.23 that we obtain (assuming / = 20°) is con- 
sistent with the value of 0.22 ± 0.06 obtained byl Miller et alJ 
(|2006), and this may indicate little change in the system ge- 
ometry even though we are observing at a lower flux. How- 
ever, these values of ft/lTT are much less than the values of 
unity or even larger that have been seen for G X 339-4 as well 
as other black hole systems in the SPL state dZdziarski et al.l 
1^03; Miller etal. 2004). Although we are including rel- 
ativistic effects that were not included in earlier work by 
IZdziarski et alj ( 120031) . our covering fraction and power-law 
index (F) values are very similar to the values IZdziarski et alj 
(l2003h used to demonstrate a systematic drop in the covering 
fraction as the source hardens. Thus, these results are consis- 
tent with a significant change in geometry between the SPL 
and the hard state, but not in the hard state between 0.056 and 
0.008 Ledd- In addition, we note that geometry is not neces- 
sarily the only change that occurs between the SPL and the 
hard state. 

The inner disk radii implied by the reflection model param- 
eters a lso suggest l i ttle or no change in ri„ down to 0.008 Lndd- 
While lMiller et all (|2006|) obtain nn = 4.0 ± 0.5, our reflection 



fits indicate maximum (90% confidence) rin values of 5.0 for 
for Spectra #1 and #2, and the spectra are consistent with lit- 
tle change from r^^ ^ 4 over the 0.008-0.056 Lsdd luminosity 
range. However, for Spectrum #2, since the iron line is not re- 
quired at high statistical confidence, we re-fitted the spectrum 
without the emission line, allowing only the blurred reflection 
component to constrain rin. In this case, we obtain = 4. 
so that we should not rule out the possibility that the inner 
disk radius increases to R^ at the lowest luminosity that 
we are sampling. 

Of course, these radius constraints are only valid if the blur- 
ring that is clearly present is caused by the Doppler boost- 
ing due to the motions of the material around the black 
hole and the gravitational redshift, and it is worthwhile to 
consider whether other physical effects could cause similar 
blurring. Probably the best developed competing model at- 
tributes the reflection component to Compton downscattering 
of the so urce's X-ray emission in a large-scale and power- 
ful wind (iTitarchuk & Shrade3l2005|) . The model requires 
hard X-ray emission with a power-law photon index of F < 2, 
and the model has been shown to be able to reproduce the 
sh ape of the ~20-200 keV ha rd state spectra from GX 339- 
4 (ITitarchuk & Shrade3l2005h . More recently, it was shown 
that this Comptonization model can produce blurred and red- 
shifted iron lines due to fluorescence of the wind material 
(|Laurent & Titarchuk 2007). Although this full model is not 
currently available for fitting our GX 339^ spectra, in the 
future, it would be interesting to see if it can explain all the 
spectral features we observe. 

5.2. The Thermal Component in the Hard State 

Further evidence for an optically thick disk in the hard state 
that is not highly truncated comes from the presence of a ther- 
mal component that has now been seen in the hard state spec- 
tra of several black hole systems. In our case. Spectrum #1 has 
a significant thermal component while the thermal component 
may be present in Spectrum #2, but it is not required at very 
high significance. Thus, we focus on comparing our Spec- 
tru m #1 param e ters to the disk-blackbody parameters found 
bv lMilleretalJ (12006) . As shown in Tabled the inner disk 
temperature of 0. 193 ± 0.012 keV for Spectrum #1 is signif- 
icantly lower than the value of kT^^ = 0.39 ± 0.04 keV mea- 
sured when the 3-100 keV flux of the source was 2.4 times 
brighter While a drop in temperature does not necessarily 
signal a change in accretion geometry, the disk-blackbody 
normalization is related to the disk inner radius according to 
Mjbb tx so the fact that A^dbb is significantly larger for 
Spectrum #1 (see TableO could indicate an increase in R^ji. 

To further investigate the difference between the values of 
Mjbb in the two spectra, we note that for these low temper- 
atures, the disk component is strongly impacted by interstel- 
lar absorption, and the value of A^h that we obtain i s some - 
what higher than the value derived by iMiller et"a]l (l2006l) . 
Thus, we refitted Spectrum #1 after fixing the column den- 
sity to the ' Miller et al.l (l2006h value. We obtain only a slightly 
higher temperature of kTi^ = 0.220 ± 0.013 keV and A^dbb = 
6100^[50Q. Thus, A^DBB is still nearly an order of magnitude 
higher fo r Spectrum #1 com pared to the value of ^700 ob- 
tained by lMiller et alJ (l2006l) . Another (not independent) way 
to examine the question of whether R^ changes is that the 
disk-blackbo dy flux (Fdisk) shou ld be proportional to kT^ for 
constant R^ (iMitsuda et alJll98 4). Based on the parameters 
with A^H fixed, the bolometric disk-blackbody flux for Spec- 
trum #1 is 3.2 X 10"'" ergs cm"^ s"', while the flux from the 
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iMiller et af] (l2006h parameters is 3.4 x 10"'° ergs cm"^ s~\ 
While the fact that these fluxes are nearly the same despite a 
drop in ^Tjn from 0.39 keV to 0.22 keV could be explained 
by a change in Ri^, other physical effects could also be im- 
portant such as possible changes in the "color correction fac- 
tor" (I Shimura & T akahara 1995) or changes in the tempera- 
ture profile in the disk. It is also worth noting that the ob- 
servations were made with different soft X-ray instruments 
(XMM-Newton vs. Swiji) and differences in calibration could 
be important. However, despite these other possibilities, if the 
correct explanation is a change in the inner radius, R\a would 
need to change by a factor of ^3. 

Our results for the evolution of the disk-blackbody pa- 
rameters a ppear to b e in contrast to results obtained re- 
cently by Rvkof f et alj ([2007D for another accreting black hole 
(XTE J1817-330) in the hard state. Using -20 Swift obser- 
vations covering inn er disk temperatures from 0.2 to 0.8 keV, 
iRvkoff et al.l (l2007h found a relationship close to Fdisk oc kT^ 
(they actually found an index of 3.3 ± 0.1 or 4.3 ±0.1 de- 
pending on the model they used for the non-thermal com- 
ponent). Although XTE J1817-330 does show some de- 
viations fr om the trend for in dividual data points (see Fig- 
ure 3 from IRvkoff eTaH llool . none of the deviations from 
the Fdisk oc kT^ trend are as large as we see for GX 339^, 
suggesting that we may be seeing a different evolution for 
GX 339-4. 

5.3. Implications for the Hard State Geometry 

Our spectra of GX 339-4 provide evidence for an optically 
thick accretion disk in the hard state. At L/Lem = 0.023, the 
evidence comes from both blurred reflection and iron line fea- 
tures as well as a significant thermal component. At L/Lndd = 
0.008, the evidence primarily comes from the blurred reflec- 
tion component. The most constraining measurement of r^ 
comes from the reflection component, and these measure- 
ments require < 5 at 0.023 Lndd and ri„ < 10 at 0.008 ^Edd- 
If the value of ri,, = 2-3 measured in the SPL state represents 
the location of the ISCO, then our hard state measurements 
imply a change of no more than a factor of 2.5 and 5 (for the 
two luminosities, respectively) greater than the radius of the 
ISCO (although see caveats discussed above). 

Recent theoretical work suggests a possible geometry that 
may be consistent with these observations. It is found that 
two physical processes can be important to causing material 
to condense out of an ADAF, leaving an inner optically thick 
disk. First, as an ADAF is forming, the material close to the 
ADAF/disk boundary will be significantly cooled via conduc- 
tive cooling, and w ill cause ADAF material to recondens e 
back into the disk dMever. Liu & Mever-HofmeisteJl2007h . 
Secondly, the soft photons from the optically thick disk can 
Compton-cool the ADAF, which also leads to condensation 
(iLiu et al .1120071) . For a relatively large range of mass accre- 
tion rates and viscosity parameters, these effects lead to inner 
and outer optically thick disk s with an ADAF filling the gap 
in between. As discussed in ILiu et al.l (12007 ). such a model 
would apply for the brighter portion of the hard state and 
could explain the small inner radii inferred from reflection fits 
in the hard state and the ^^0.2-0.3 keV thermal components. 
This geometry is also at least qualitatively consistent with the 
observed H/Itt values. 

Although such condensation can occur for relatively brigh t 
portions of the hard state, the calculations of Liu et aP (l2007h 
still indicate that below L/Lem ~ 0.001, the inner disk will 
evaporate, leaving only ADAF inside some truncation radius. 



Thus, it is notable that our observations of GX 339^ as well 
as most, if not all, of the cases where small inner disk radii 
have been inferred for the hard state via reflection modeling 
or the presence of soft components hav e occurred above thi s 
level (Miller. H oman & MiniuttH 120061: 'Rykoffet al.' '2007). 
At lower luminosity, there have been observations that do 
not necessarily require soft components in the X-ray band. 
For example, for the b lack hole system XTE J1118H-480, 
iMcChntock et al.l (1200 Ih find evidence for a soft component 
at the very low temperature of —24 eV when the system was 
near 0.001 Ledd (assuming a dis tance of 1.8 kpc and a black 
hole m ass of 7M0). In addition. iTomsick. Kalemci & Kaarel 
(l200l observed XTE J 1650-500 at levels of lO'^'Ledd (as- 
suming a distance of 4 kpc and a black hole mass of IOMq), 
and did not detect a soft component or iron features. Finally, 
very high quality spectra have been obtained for several sys- 
tems in quiescence L/Lem ^ lO"*" or lower wit hout evidence 
for a thermal component or iron featur es (e.g.. lBradlev et all 
l2007HCorbel. Tomsick & Kaaretll2006h . 

5.4. Implications for the Compact Jet 

The possibility of an inner optically thick disk in the hard 
state has very interesting implications for the production of 
compact jets in the hard state. For GX 339-4, we detect the 
compact jet in the radio band contemporaneous ly with the 
times that we obtained Spectra #1 and #2, and Miller et al.l 
(l2006l) also report the presence of a compact jet during their 
hard state observation. Cygnus X-1 provides another ex- 
ample of a bright hard state black hole with a compact jet. 
These examples imply that a compact jet can be produced 
when the inner optically thick disk is present. Furthermore, 
it is notable that fitting multi-wavelength spectral energy dis- 
tributions from radio-to-X-ray observations of GX 339^, 
Cygnus X-1, and another black hole system, GRO J 1655^0, 
with a compact jet model give values of 3.5-10 for the 
radius of the jet at its base dMarkoff. Nowak & Wihnslllool 
iMigliari et al.l2007h . These small radii are consistent with the 
jet being launched at or within the inner edge of the disk, and 
may imply that the production of the compact jet is closely 
linked to the inner disk. 

However, at the same time, radio observations show that 
compact jets can also be pro duced at very low luminosi- 
ties by quie scent black holes (iGallo. Fender & Hvnesll2005t 
iGallo et an F2006). This brings the role of the inner disk in jet 
production into question since it is unclear whether the inner 
disk persists to these low luminosities. As discussed above, 
the conductive and C ompton cooling mechanisms considered 
by iLiu et al.l (l2007h indicate that the disk should evaporate 
below 0.001 Lfidd- However, it may still be worth consider- 
ing whether other cooling mechanisms can maintain the inner 
disk to lower levels. For example, Meier (2005) has explored 
the possibility of Magnetically Dominated Accretion Flows 
(MDAFs), and there are indications for an inner disk region in 
the hard state where magnetic cooling is important. Clearly, 
more theoretical studies as well as observations at these very 
low luminosities are important for a full understanding the full 
set of conditions that are required for compact jet production. 
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TABLE 1 
Observations of GX 339-4 



Satellite 


Observation 


Date in 


Start Time 


Stop Time 


Exposure" 




ID (ObsID) 


2007 


(UThour) 


(UT hour) 


(s) 


Spectrum #1 


Swift 


00030943001 


May 25 


17.10 


23.80 


6150 


RXTE 


92704-04-02-00 


May 25 


17.99 


20.13 


4976 


RXTE 


92704-04-02-02 


May 25 


21.14 


21.69 


1808 


Spectrum #2 


Swift 


00030943002 


June 10 


16.70 


21.70 


4685 


Swift 


00030943003 


June 12 


1.29 


21.98 


5509 


Swift 


00030943004 


June 14 


1.01 


22.38 


3204 


RXTE 


92704-03-28-00 


Jime 10 


4.43 


5.63 


2320 


RXTE 


92704-03-29-00 


June 1 1 


8.44 


9.51 


3392 


RXTE 


94704-03-29-01 


June 12 


3.27 


3.64 


1056 


RXTE 


94704-03-30-00 


June 13 


6.08 


6.64 


1952 


RXTE 


94704-03-31-00 


June 14 


15.29 


16.28 


2112 



"This is the exposure time on the target (GX 339^) obtained by the XRT instrument (for Swift) or the PCA instrument (for RXTE). 
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TABLE 2 

GX 339-4 Spectral Fits with Basic Models 



Model" 


(10-' cm-2) 


r 




(keV) 






Spectrum #1 


PL 


3.1±0.1 


1.66 ±0.01 


0.113 ±0.002 






651.8/213 


PL+DBB 


8.5 ±0.7 


1.69 ±0.01 


0.125 ±0.003 


0.178 ±0.007 


76,000_25ooo 


532.0/211 


SMxPL 


2.6±0.1 


1.58±0.01 


0.104 ±0.002 






444.7/211 


SMx (PL+DBB) 


7.2 ±0.7 


1.59 ±0.02 


0.109 ±0.003 


203+'' '"2 


23 000+'^-™ 


338.3/209 


Spectrum #2 


PL 


4.9 ±0.2 


1.61 ±0.02 


0.0329 ±0.0008 






321.6/213 


PL+DBB 


9.0±1.1 


1.64 ±0.02 


0.0367 ±0.0013 


0.157±0.011 


35 OOO^'"-™ 


280.1/211 


SMx PL 


4.2 ±0.2 


1.50 ±0.02 


0.0294 ±0.0008 






210.6/211 


SMx (PL+DBB) 


7.0±1.4 


1.52 ±0.03 


0.0315±0.0015 


0.184:«;«i 


5 900+* ™" 


200.2/209 



fits, we fixed the width of the smeared iron edge to 10 keV. 
''En'ors on all parameter values are 90% confidence (Ax^ = 2.7). 
'^Units are photons cm -2 S-' keV-' at 1 keV. 



TABLE 3 

GX 339-4 Spectral Fits with Reflection 



Parameter" Spectrum #1 Spectrum #2 2004 Spectrum* 



(10-' cm"^) 


7 1+"'' 
'■^-0.8 


7 3+' -^ 
'■-^-1.3 


3.7 ±0.4 


r 




'■°-'-0.03 


1.41 ±0.03 


A'pl" 


n 1 1 0+0.008 

"■'^'^-0.006 




0.32 ±0.03 


«:7i„ (keV) 


0.193±0.012 


0.161+0017 
14 000+^^'^^° 


0.39 ±0.04 


A'dbb 


29 000+^' "™ 


700 ± 200 


i (degrees) 


20.0 


20.0 


20+^ 


Vi/2iz^ (pexriv) 


A 99+0.06 
'^•^^-0.05 


24+<'" 


0.22 ±0.06 


(ergs cm"' s"') (pexriv) 


10 000+^ °™ 

1U,UUU 5 


7 000+**'°™ 
' ' """-4.000 


1,000 


£line (keV) 


6 9+'' 2 


6 7+0 '' 

"■'-0.3 


6.8±0.1 


A'line (photons cm"^ s"') 


(7.4 ± 4.4) X lO"'' 


(2.4+1^) X 10"^ 


(3.5 ± 0.3) X 10"^ 


Rin/Rg (kdblur) 


3 6+' '' 


2 9+2' 


4.0 ±0.5 


/(kdblur) 


■X 9+0.5 
-'■^-0.6 


3.1 ±0.4 


3.0 




324.7/205 


191.3/205 


2120.5/1160 


Absorbed Flux (ergs cm"- s"')^ 


2.1 X 10"' 


7.4 X 10"'0 


5.4 X 10"' 


Unabs. Flux (ergs cm"^ s"')'' 


2.2 X 10"' 


7.7 X 10"'0 


5.5 X 10"' 


Luminosity (ergs s"' )' 


1.7 X 10" 


5.9 X 10'* 


4.2 X 10" 


Luminosity (LEj^y 


0.023 


0.008 


0.056 


Iron Line Equivalent Width (eV) 


140 ±90 


■ 40+150 


~160 



"Errors on all parameter values are 90% confidence (Ax^ = 2.7). 

*These are the parameters obtained bv lMiller et al] (2003) for the 2004 XMM-Newton and RXTE observations of GX 339^. 
^Units are photons cm"- s" keV"' at 1 keV. 
''Reflection covering factor. 
'Ionization parameter. 

^The power-law index for the radial emissivity profile. 
•"Absorbed flux in the 1-100 keV band. 
''Unabsorbed flux in the 1-100 keV band. 
' 1-100 keV luminosity assuming a distance of 8 kpc. 

^1-100 keV luminosity in Eddington units assuming a distance of 8 kpc and a black hole mass of 5.8Mq. 



